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ABSTRACT 

We present a first look at the local luminous infrared galaxy IRAS 04296+2923. This barred spiral 
galaxy, overlooked because of its location behind the Taurus molecular cloud, is among the half dozen 
closest (D= 29 Mpc) LIRGs. More IR-luminous than either M82 or the Antennae, it may be the best 
local example of a nuclear starburst caused by bar- mediated secular evolution. We present Palomar 
J and Pa j3 images, VLA continuum maps from A =20-1.3 cm, a subarcsecond Keck LWS image at 
11.7 jam and Owens Valley Millimeter Array CO(1-0), 13 CO(1-0), and 2.7 mm continuum images. The 
J-band image reveals a symmetric barred spiral galaxy. Two bright, compact mid-infrared and radio 
sources in the nucleus mark a starburst that is energetically equivalent to ~ 10 5 07 stars, separated 
by ^ 50 pc. This is probably a pair of young super star clusters, with estimated stellar masses of 
~ 10 T Mq each. The nuclear starburst is forming stars at the rate of ^ 12 + 6 M yr _1 , or about half 
of the total star formation rate for the galaxy of ~ 25 ± 10 M yr _1 . IRAS 04296+2923 is very bright 
in CO, and among the most gas-rich galaxies in the local universe. The 12 CO luminosity of the inner 
half kpc is equivalent to that of the entire Milky Way. While the most intense CO emission is extended 
over a 15" (2 kpc) diameter region, the nuclear starburst is confined to within 1-2" (150-250 pc) of 
the dynamical center. Based on masses obtained with 13 CO, we find that the CO conversion factor 
in the nucleus is higher than the Galactic value, X^q by a factor of 3-4, typical of gas-rich spiral 
nuclei. The nuclear star formation efficiency is nuc M gas /SFR nuc = 2.7 x 10 -8 yr _1 , corresponding to 
a gas consumption timescale, r^p c ~4x 10 T yrs. The star formation efficiency is ten times lower in 
the disk, with rj^ fe ~ 3.3 x 10 8 yrs. The low absolute star formation efficiency in the disk implies 
that the molecular gas is not completely consumed before it drifts into the nucleus, and is capable of 
fueling a sustained nuclear starburst. IRAS 04296+2923 appears to be beginning a 100 Myr period 
as a LIRG, during which it will turn much of its 6 x 10 9 M of molecular gas into a nuclear cluster 
of stars. 

Subject headings: galaxies individual (IRAS 04296+2923,2MASX J04324860+2929578) — galax- 
iesistarburst — galaxies: ISM — radio lines — radio continuum: galaxies 



1. INTRODUCTION 



Luminous (L/# > 10 11 L ) and ultraluminous 
(Lir > 10 12 L ) infrared galaxie s are powered by prodi- 
gious amounts of star formation ([Joseph fc Wri ght 1985; 
Con don et al. 1991 ; Sanders & Mirabel 1996 ;|GenzeI et ah] 
1998^ Lutz et al.|1998||Downes Solomon||1998|. While 
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U LIRGs ar e typically associated with mergers ( |Sanders 
|et al.|[l988| ), in LIRGs, star formation may be generated 
by interactions or secular evolu tion, often related to bars 
( Kormendy 



Ken nicutt 
(poller et al.|1987p 



The" 



2004 ) . LIRGs are relatively rare 
:ew LIRGs close enough for their 



internal structure to be resolved are valuable targets for 
the study of the local physics of the starburst /gas inter- 
action, such as when and where the star formation takes 
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place within a galaxy. 

IRAS 04296+2923 (hereafter IRAS 04296) lies behind 
Ay ~ 5.5 mag of extinction from the dark cloud LI 500 
in Taurus (Table pi). Low-resolution optical and HI spec- 
tra allowed it to oe identified as a galaxy (IStrauss et al.l 



1992||Chamaraux et al 1995| ). IRAS 04296 has an IK. 
luminosity of 9.8 x 10 ltr T^(29 Mpc, for H 



_1 Mpc~ 



71 km 

Sanders et al.|2003 ), making it the fifth most 



IR-luminous galaxy within 30 Mpc, at the lower bound 
of the LIRG class. It is 60% more IR-luminous than 
M82, and 35% more than the Antennae, exceeded only 
by NGC 1068, NGC 2146, NGC 7552, and NGC 1365 in 
the local universe. 

Little is known about IRAS 04296. We discovered its 
nuclear starburst in a subarcsecond mid-infrared imaging 
survey with the Long Wavelength Spectrometer (LWS) 
on the Keck Telescope; it has a warm IRAS 60/100/im 
color that is a diagnostic of concentrated star form ation . 
The bright double mid-IR source revealed by LWS (j |3.1| ), 
< 1" (150 pc) across, suggested the presence of a new 
starburst galaxy in the local universe. 

In this paper we present near-infrared (NIR), mid in- 
frared (MIR), millimeter and centimeter radio continuum 
images of IRAS 04296, as well as NIR and millimeter 
spectral imaging. While its location in Taurus hinders 
optical and ultraviolet observations of IRAS 04296, in- 
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frared and radio observations are unaffected, and they 
are the most reliable probes of star formation in dusty 
starbursts. 

2. OBSERVATIONS 

2.1. Near- Infrared J -Band and Pa f3 Images from the 
Hale Telescope 

IRAS 04296 was observed using the Palomar Hale 5- 
meter telescop on 2005 Januar y 23 UT with the W ide- 
field Infrared Camera (WIRC; | Wilson et aT] [2003b at 
prime focus, in broad-band J and narrow- band Pap fil- 
ters. The camera has a pixel size of 0'/25. Exposure 
times on-source were 11 minutes (3 x 20 sec frames, 11 
pointings) for the continuum band and 18 minutes (3 
x 60 sec frames, 6 pointings) in the emission-line band, 
with equivalent off-source exposure times to measure the 
sky. The observing conditions were good with seeing at 
l'/O FWHM. Median dark current and filtered sky images 
created for each band were subtracted from the frames. 
Finally a coadded mosaic was produced from each set of 
on-source frames. We scaled the counts in each band's 
mosaic, via aperture photometry of several bright, un- 
crowded stars in the field around the galaxy, and sub- 
tracted the J-band mosaic (see Fig.JlJi) from the Pa/3 mo- 
saic, to produce an approximate continuum subtracted 
image (see Fig. [!}>). The images in the broad-band and 
narrow band filters were not explicitly flux-calibrated us- 
ing observations of standard stars, so we are unable to 
derive a value for the net Pa/3 flux in the continuum sub- 
tracted image. 

2.2. Mid-Infrared Images from Keck 

IRAS 04296 was observed on 24 December 2002 with 
the Long Wavelength Spectrometer (LWS) in imag- 
ing mode on the Keck I Telescope, with the standard 
11.7/im/SIC filter. The standard stars were a Aur, 
a Aries and j3 Leo, observed both before and after the 
source exposures. High quality images at both wave- 
lengths were obtained (Fig. J2|; however, the absolute 
flux calibration was uncertain due to intermittent in- 
strumental problems, and so we reobserved the galaxy 
on 31 January, with the same parameters as confirma- 
tion. Absolute flux calibration is estimated to be good 
to 10%, based on derived flux counts on standard star 
before/after observations throughout the night. Stellar 
images were about 5 pixels or 0'/4 FWHM. 

2.3. Radio Continuum Images from the VLA 

The radio data were taken with the NRAO VLA0 
There were two projects: AB1092 on 2003 July 7 with the 
A configuration for high angular resolution at A = 20 cm, 
6 cm, 3.5 cm, 2 cm, and 1.3 cm (L,C,X U and K bands) 
and AT309 on 2005 July 5 with the array in BC con- 
figuration for moderate resolution at 1.3 cm only. 3C48 
was the primary calibrator for both runs. The 20, 6, and 
3.6 cm data were processed with AIPS, including baseline 

7 Based on observations obtained at the Hale Telescope, Palo- 
mar Observatory, as a part a continuing collaboration between the 
California Institute of Technology, NASA/JPL, and Cornell Uni- 
versity. 

8 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



corrections. The 2 cm and 1.3 cm data were observed in 
fast-switching mode and reduced according to the AIPS 
prescriptions for high frequency calibratio n. T he details 
of the radio imaging are discussed below ({3.2) 



2.4. Millimeter CO Lines and Continuum from OVRO 

Simultaneous aperture synthesis observations of the 
12 CO(1-0) transition (115.271 GHz) and the 13 CO(1-0) 
transition (110.201 GHz) of IRAS 04296 were made with 
the Owens Valley Radio Observatory (OVRO) Millimeter 
Interferometer between 2003 November 1 and 2004 May 
8. The interferometer consisted of six 10.4 m antennas 



with cryogenically cooled SIS receivers (Scoville et al. 
|1994[ ). Observations in "Compact", "Low" and "High" 
configuration were obtained, with system temperatures 
(single sideband) ranging from 220 - 450 K at 115 GHz. 
64x4 MHz channels were used to cover both transitions, 
giving a velocity resolution of 10.5 km s _1 for CO(1-0). 
The phase center is within 0'/5 of the 2 /am peak listed 
in Table [l] and the adopted systemic velocity (LSR) was 
2100 km s -1 . Simultaneous low spectral resolution, wide- 
band (128x31.25 MHz) COBRA observations were used 
to generate a 3 mm continuum image (y Q ~ 111 GHz), 
free from 12 CO(1-0), CN(l-O) and 13 CO(1-0) contamina- 
tion. The data were calibrated using the MMA software 
package. Phase calibration was done by observing the 
quasar J0336+323 every 25 minutes. Absolute flux cali- 
bration was done using Uranus as primary flux calibrator 
and 3C273 and 3c454.3 as secondary flux calibration, and 
is good to ~10%. 

Mapmaking was done in MIRIAD and subsequent data 
analysis and manipulation was done with the NRAO 
AIPS package. Maps were generated with natural 
weighting to maximize sensitivity, uniform weighting to 
maximize resolution and robust = weighting to simul- 
taneously optimize sensitivity and resolution. Unless 
stated the analysis results from the naturally weighted 
data. Integrated intensity images are straight moment 
maps with all emission brighter than 1.3 a (la for 13 CO) 

per channel included. The OVRO primary beam is ^64 
at 115 GHz. Corrections for the primary beam have not 
been applied so emission at the edge of the primary beam 
will be somewhat underestimated. Structures extended 
on scales larger than ~52 for CO (1-0) will be resolved 
out by the interferometer. No single-dish observations 
of this galaxy exists, so no estimate of the amount of 
resolved-out flux is possible for the mm- wave images. 
However, 52 corresponds to spatial scales of 7.3 kpc 
and given that the emission must be extended on this 
scale in an individual channel to be missed, we consider 
it unlikely that significant flux is resolved out. If large 
amounts of flux are missed then the observed molecular 
richness is a lower limit. 

3. RESULTS: A FIRST LOOK AT IRAS 04296+2923 

Near-infrared images of IRAS 04296 reveal a normal, 
barred spiral galaxy in starlight. The bright, concen- 
trated nuclear starburst is prom inent at MIR and radio 
continuum wavelengths. In §3.1| we discuss the large field 
Palomar near-IR images of the galaxy and its Pa/3 emis- 
sio n, an d the Keck high resolution 11.7/im MIR image. 
In £3.2 we discuss the VLA radio images, which were 
mapped in four ways, one to match beams for extended 
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emission dominated by synchrotron emission, and two for 
matching compac t em ission dominated by free-free from 
HII regions. In §3.3| we discuss the CO and molecular 
gas. 

3.1. Infrared Views 

The Palomar J band continuum image and a narrow 
band continuum-subtracted image of Pa/3 in IRAS 04296 
are shown in Figure [T] The near-infrared continuum at 
J-band reveals the stellar population. IRAS 04296 is an 
inclined, barred spiral with outer ring-like spiral arms 
(SBb/c(r)). The total extent of the emission is roughly 
2-2.5', or 17-20 kpc (8-10 kpc in radius). The stellar 
distribution appears bi-symmetric and undistorted. The 
Paj3 image shows line emission to be very concentrated 
and bright in the galactic nucleus, with some very weak 
emission along the inner arms. The compact nature 
of the bright Pa/3 emission shows that the starburst is 
highly concentrated at the nucleus. Since the images 
are not photometric, and extinction is probably variable 
across the galaxy, we cannot precisely quantify the rela- 
tive brightnesses of disk and nuclear star formation from 
the Pa/3 image. 

Mid-IR continuum emission traces warm (~ few hun- 
dred K) dust heated by young stars in the starburst. 
The high resolution (0.3") mid-IR continuum image at 
11.7/xm made with LWS on Keck of IRAS 04296 is shown 
in Figure [5] the field of view of this image is only 10" x 
10", or about 1.4 kpc on a side. A similar 18.75/im im- 
age (not shown) will be published in a subsequent paper 
(Turner et al. 2010, in prep.). The 11.7/xm filter contains 
two PAH features; however, based on observations of the 



emission from HII regions is compact with sp ectra that 



are eithe r flat or rising with frequency (e.g., [Turner fc 
Ho|1994[ ). Except at 1.3 cm, the robustly- weighted VLA 



starburst galaxies NGC 7 714 and Arp 220 (Brandl et al. 
2004; |Armus et aT1|2007|) we estimate that at most 15% 
ol the 12/xm continuum is due to PAH emission. 

The mid-IR emission comprises a bright source, which 
cannot be resolved at the 0"3 diffraction limit of Keck 
at 11.7/im, and a secondary source or tail extending to 
the east of the main source, separated from the main 
source by ~1"5 (~200 pc). The total flux density is 
^ii.7/im = 680 mJy. These two compact, mid-infrared 
sources account for roughly half of the IRAS 12 jira flux 
of the entire galaxy. 

3.2. Radio Views 

VLA maps of IRAS 04296 across the cm waveband are 
presented in Figure [3j with fluxes and noise values given 
in Table [U The radio continuum emission is a combina- 
tion of thermal bremstrahlung emission from the HII re- 
gions and nonthermal synchrotron emission from super- 
novae and their remnants associated with the starburst. 
The brightest radio emission, like the mid-IR emission, 
is concentrated within the central 1-2" region of ~ 150- 
250 pc extent; higher resolution maps shown below give 
a better comparison with the Keck mid-IR images. 

Separating free-free emission from the HII regions 
from the synchrotron emission is possible because these 
sources have different spectral and spatial characteristics. 
However, complications arise when comparing aperture 
synthesis maps at different frequencies, so the spectral 
analysis must be done with care. The synchrotron com- 
ponent in galaxies is spatially extended with a spectrum 
that falls steeply with increasing frequency. Free-free 



maps oi Figure [3] were done with a single array configu- 
ration, so the maps have responses to extended structure 
that vary with frequency; shorter wavelength maps are 
less sensitive to extended emission. The shortest base- 
lines determine the maximum detectable spatial scales, 
@max ~ ^/Bmin, which are listed in Table [l] So the 
maps of Figure [3] cannot be compared directly to obtain 
spectral index maps. Since single dish or lower resolu- 
tion VLA fluxes do not exist for IRAS 04296 except at 
20 cm, we cannot estimate the effects of missing short 
baselines except at 20 cm. Our A array maps, with total 
flux S20 cm = 140 zb 7 mJy, recov er the entire 20 cm flux 
recorded by |Condon et al | ( |1996[ ) , and so the 20 cm map 
is probably a good representation of the total radio emis- 
sion, although there may be faint extended emission that 
is undersampled. From the 20 cm image we can see that 
the extended radio emis sion is aligned along the apparent 
bar of the galaxy (£3.3). The bright mid-IR sources and 
the peak of the radio emission mark a starburst located 
a t the inner end of the northwestern arm traced by CO 

Free-free and nonthermal synchrotron sources can be 
separated by determining the spectral index, a, of the 
emission (S oc v a ) so long as the beams are matched, 
not only in beamsize (6 m i n ), but also in 6 ma x, defined by 
Umin = B^ n j\ v min B^fjX. Toward this end we 
have made a series of "cut" maps, which are high-pass 
spatially filtered images with common min and max ^ 
that enable us to construct images with nearly identi- 
cal (u, v) coverages and allow us to compare the radio 
spectra for the compact sources, at the cost of losing 
response to extended emission. Since the star formation 
arises largely in compact sources <1" (150 pc), we expect 
to detect nearly all of the free-free emission in the "cut" 
maps and resolve out most of the synchrotron emission, 
which is our goal. 

The first set of "cut" maps were made to match the 
beams at 20, 6, 3.6, and 1.3 cm, with a consistent 
(u,v) m i n of 18k\ (6max ~ 10", 1.4 kpc). The beam 
for this series of maps is 0"82 x 0"29, p.a. -89°. For 
the 1.3 cm emission, we used the BC configuration data. 
From the cut maps, we mapped the spectral index a, 
shown in Figure [4] The spectral index is negative over 
the entire range, confirming that the extended emission 
cm- wave emission is largely no n-thermal, typical of large 
spiral galaxies and LIRGs ( e -ff- Condon et al.|199l] |Con-| 



|don||1992| |Clemens et al . 

maps are measured tor a 



2008 ). Fluxes Irom this seToT 
3" box (630 pc 



420 pc) 

are shown in Table [TJ 

A second set of "cut" maps were made from the shorter 
wavelength data to image the compact emission at 3.6, 2 
and 1.3 cm. For these maps, (u,v) m i n = 50 kA (0 max ~ 
4", or ~ 550 pc) and they were convolved to match the 
/ /27x0 / /23, p.a. -90° beam (0 min ) of the 3.6 cm map. 
These cut maps are shown in Figure |5j and fluxes are 
given in Table [T] The peak flux densities are the same 
within observational uncertainties for all three maps, 1.1- 
1.5 mJy/beam. The (u,v) restriction has eliminated 
about two-thirds of the total 3.6 cm flux as compared 
to the map of Figure [3| The compact emission that re- 
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mains in these maps has a nearly flat spectrum: this is 
emission from HII regions. These radio images are sim- 
ilar to the 11.7/im image of Figure [2j given the lower 
resolution (0'/3) of the mid-IR image. We identify the 
brightest, elongated radio source with the mid-IR peak, 
and the extension visible to the southeast in the 3.6 cm 
image (Fig. [5]) with the extension in the mid-IR image. 

A third set of "cut" maps, designed to match beams at 
2 and 1.3 cm, give the highest resolution images. Like the 
previous set of maps, they are cut to (u,v) m i n = 50 kA 
but these maps are instead convolved to match the 
smaller 2 cm beam of 0'/ 15 x0V 14, p. a. -66°. These 
images resolve the brightest continuum source (Fig. [6| 
into a double source separated by 0'/2, or 30 pc oriented 
east-west. The eastern source is the stronger. Slight 
differences in the spatial structure of the double source 
at the two frequencies might be due to differences in az- 
imuthal (u,v) coverage or to 1.3 cm seeing. This source is 
flat in radio spectral index, indicating that the compact 
emission is largely free- free emission from HII regions. 
Comparison of the three sets of maps suggests that the 
total amount of flat-spectrum, compact emission in the 
nucleus of IRAS 04296 is -10 mJy. 

Continuum emission at 2.7 mm was detected with 
OVRO (Figure [7]i). The peak 2.7 mm continuum inten- 
sity is 10±1 mJy bm" 1 (4. 6x3. 7) with a total detected 
flux of 11±2 mJy (Table IJ. The millimeter continuum 
emission is confined to the nuclear starburst, with po- 
sition and extent consistent with the (uncut) cm-wave 
and infrared continuum. At this wavelength, potential 
continuum sources include synchrotron, free-free emis- 
sion and dust. Extrapolation of synchrotron emission all 
the way from 20 cm to 2.7 mm is dangerous since it is 
sensitive to the exact spectral index, but using the 20 cm 
emission seen in the first 'cut' map and a =-0.8 implies 
a contribution o f 2-3 mJy. Based on our estimated nu- 
clear gas mass (j |4.2| ) we expect that —1-2 mJy of the 3 
mm continuum is from dust emission. From this, we con- 
clude that the millimeter continuum flux associated with 
star formation is >5-7 mJy at 2.7 mm, consistent with 
the fraction of the flux from the compact flat spectrum 
sources seen directly in cm-wave continuum. 

3.3. Molecular Views 

The first observations of molecular gas towards IRAS 
04296 are shown in Figures [7 48] The CO (1-0) integrated 
intensity (Fig. [7]) is remarkaMy bright, extending over a 
2.5 kpc (diameter) region, peaking at the nucleus. (When 
no superscript appears the most abundant from, 12 C 16 
is implied; ' 13 CO' refers to the less abundant 13 C substi- 
tuted isotopologue.) The CO emission extends well be- 
yond the Pa f3 emission region (Figure [lb) and the radio 
continuum (Figure [8]). Beyond the nucleus the molecu- 
lar gas is bar-like, extended southeast-northwest. Weak 
emission is also seen to the north and west of the nucleus. 
The CO morphology closely follows the stellar morphol- 
ogy seen in the Palomar J band and 2MASS K s images, 
as shown in Figure [7] 

Towards the nucleus, antenna temperatures in the ro- 
bustly weighted image (Figure (8| peak at 10 K averaged 
over the inner 160 pc radius. This is a remarkably high 
brightness temperature for a galaxy at the distance of 
IRAS 04296; it implies that the beam filling factor for 



the CO must be near unity over regions 300-400 pc in 
diameter. The CO line widths are broadest towards the 
nucleus but are fairly modest (250 k m s -1 , FWHM) com- 
pared to other luminous starbursts fl S olomon , Downes~fc 



Radford 1992; |Solomon et al.|fl997| |Downes & Solomon 



1998), and drop to ^50 km s towards the outer molec 



ular arms. The velocity field is consistent with a rotat- 
ing disk, except for modest perturbations (<50 km s _1 ; 
see below). We have fit the axisymmetric velocity field 
with a Brandt rotation curve; the best-fit parameters are 
given in Table [T] The separate CO component north of 
the nucleus follows the same velocity field as does the 
bar and nucleus. This suggests that these clouds are a 
continuation of the galaxy's disk (the outer spiral arms). 

The molecular morphology is bar-like. The best fit for 
the kinematic major axis from the CO velocity field is at 
a position angle of 252° (to the receding axis). This is 
perpendicular to the long axis of the CO and IR morphol- 
ogy. Therefore the observed bar-like structure seen in CO 
is not a result of inclination, but is intrinsic. The strong 
bar deduced for the stellar distribution makes it impos- 
sible to separately constrain the bar strength and the 
galaxy's inclination angle, z, over the inner arcminute. 
Instead we estimate the inclination angle from the axial 
ratio of the outer, presumably more circular, spiral arms 
seen in the J band image (Fig. [T^i). The measurement 
yields i = 50 ± 4° and is consistent with the same p. a. 
determined from the CO velocity field. Such an inclina- 
tion angle implies a peak velocity of ~190 km s _1 , or a 
dynamical mass of 1.6 x 10 9 M over the 1" (R < 500 pc) 
and 3 x 10 10 M over the radius, R = 30" (R < 4.3 kpc). 
This suggests that IRAS 04296 is similar in size to the 
Milky Way at this radius. HI observations at larger radii 
will be needed to determine a global dynamical mass for 
IRAS 04296. 

To verify that the gas distribution and kinematics are 
consistent with response to a bar potential, we generate 
a simple analytical, weak-bar model. These models treat 
the gas dissipation by adding a damping term propor- 
tional to radial velocit y to the standard stellar b arred 
orbits (e.g. |Wada|1994} (Lindblad fe Lindblad|1994p . De- 
spite the simplicity oi the model, it matches full hy- 
drodynamical simulations with surprising fidelity (e.g. 
iLindblad fc Lindb lad 1994]). The current model is that 
ol |Sakamoto et al.| ( |l999| )' except that we use an ax- 
isymmetric potential that generates a Brandt rotation 
curve. With the analytical model we can quickly search 
barred galaxy parameter space for configurations consis- 
tent with what is observed. 

Figure 9j which shows the bar models and rotation 
curve, indicates a close agreement of models with the ob- 
served morphology and kinematics. From the best mod- 
els we find that the molecular gas distribution of IRAS 
04296 is very well reproduced by a large axial ratio bar. 
We assume that the northern portion of the galaxy is 
the near side, consistent with two fairly straight trailing 
arms on the leading edge of the bar. The observed veloc- 
ity residuals have the expected signatures of inflow along 
the downstream sides of the bar arms and weaker outflow 
on the upstream (so called "spray" region) sides of the 
arms. Both modeled radial velocities and peak observed 
velocity residuals are in the range of V r ~ 10 — 40 km 
s _1 . 
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4. DISCUSSION 

4.1. The Star Formation Rate and Luminosity of the 
Starburst in IRAS 04296 

From the observed free- free emission (j |3.2| ) we can in- 
fer the Lyman continuum rate of the starburst, its Ljr, 
and star formation rate. The compact cm- wave emis- 
sion of IRAS 04296 consists of an east-west elongated 
source, centered at RA=04 /l 32 m 48!6, Dec=29°29'57". 
The brightest emission corresponds to a region O'MxO'.^ 
in size (50 pc x 30 pc) with a flux of ~3 mJy and a 
flat spectrum, located about 1" (140 pc) WSW of the 
dynamical center. There is a halo of emission, stronger 
to the south of the double source (Fig. El. The emis- 
sion measure, EM, of these HII regions is nigh; our high 
resolution maps (Figs. [5] and [6| detect only gas with 
EM > 10 6 -10 T cm -6 pc, which is characteristic of dense 
or compact HII regions in the Galaxy. The central dou- 
ble radio source has a total flux ~ 7 mJy at 1.3 cm, and 
~10 mJy at 2 and 3.6 cm. The short wavelength (2 - 0.26 
cm) continuum emission is consistent with star formation 
having a free-free flux of ~5-7 mJy, which implies an ion- 
izing flux from young stars of NL yc ~ 6 ± 2 x 10 53 sec -1 
for optically thin emission. This total Ly man contin- 
uu m rat e is equivalent to 1.4 ± 1 x 10 5 07 (Martins et 



aLl|2005t stars. Using STARB URST99 (Leitherer et al. 
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1999[ | Vazquez fc Leitherer|2005| ) and a Kroupa IMF with 
mass cutoffs ol 0.1 and 100 M , we infer a luminosity as- 
sociated with the OB stars of L OB ~ 2 ± 0.7 x 10 10 L 
for the nuclear (R<150 pc) starburst, to within a factor 
of two, given uncertainties in free- free flux and starburst 
age. The mid-IR continuum flux is consistent picture: 
the 12/im/radio flux ratio of 100 is similar to the values 
of 50-100 obser ved in both Galactic and extraga lactic 
HII regions (e.g. |Genzel et al |1982| [Ho et al.|1989| >. The 
thermal cm-wave fluxes imply a star formation rate of 
- 12 ± 6 M yr' 1 and a mass of 2 ± 1 x 10 7 M for the 
nuclear starburst. From the morphology, we infer that 
these are two massive young super star clusters, "hyper 
star clusters", each containing 50,000 O stars and with 
individual masses of ~ 10 7 M . The nuclear starburst 
contains about half of the current global star formation 
rate of IRAS 04296, based on the ratio of LWS flux for 
the nuclear source to the IRAS 12/i flux for the entire 
galaxy. 

4.2. A Super CO-Rich Galaxy: Conversion Factors and 
Estimating IRAS 04296 's H 2 Mass. 

IRAS 04296 is very bright in CO. Its total mapped CO 
luminosity is Leo — 1.4 x 10 9 K km s _1 pc 2 . The nucleus 
is the dominant feature in the CO maps. The peak inte- 
grated CO (1-0) intensity averaged over the central kpc 
(R = 3.5 ) is higher than the peak CO brightness aver- 
aged over individu al GMC scales in the nearby gas-rich 
spiral, IC 342 (e.g. [Meier Turner fc Hurtp OOO). Among 
nearby galaxies, only MGC 2146, also a L1RG, and possi- 
bly NGC 253, can rival it in the richness of its mo lecular 
emission ( jJackson fc Ho||1988| |Young et alll988[) . 

The bright CO implies a large H2 mass. For a Galactic 
conversion factor of X^ = 2 x 10 20 cm~ 2 (K km s -1 ) -1 



larger than the molecular mass of the entire Milky Way. 
This mass is comparable to the dynamical mass for the 
central kpc, as determined by the rotation curve from 
our model (Fig. [9]). Either molecular gas accounts for 
the entirety of the dynamical mass in the inner 500 pc 
radius or the CO-to-H2 conversion factor is different than 
in the Galaxy. Evidence is accumulating that the latter 
explanation is far more likely in the center of a gas-rich 
galaxy such as IRAS 04296. Studies of a variety of galax- 
ies indicate that Xf?Q as determined from the molecular 
clouds in the disk of our Galaxy is inappropriate for the 
central regions of galaxies. In (U)LIRGs it is well known 
that the Galactic conversion factor consistently overesti- 
mates (nuclear) molecular masses by factors of 3-4 (e.g. 
I Downes et al.|1993| ISolomon et al.|1997j |Yao et al.|2003| 
|Narayanan et al.||2005p . This is often called the "star- 
burst" conversion factor. However, in fact, this lower 
conversion factor also holds in gas-rich centers of normal 
spiral galaxies including relatively quiescent ones such as 



our own (e.g. Mauersberger et al.||1996[ |Dahmen et al. 
[T998l |Meier & Turner ||2001l | Weiss et al.pOOH |Pag lione 
et al.|2Q01||Meier fc Tumer||2004| |Meier, Turner fc Hurt 
2008). Since Xqo is a dynamical measure of mass, any 
contribution to the cloud linewidths from systematic mo- 
tions driven by the stellar potential will cause an over- 
estimate of the gas mass. This is particularly likely in 
galactic centers, where tidal effects on the clouds become 
important (e.g., Meier fc Turner|2004[ ) . 

Our 13 CO map provides an independent measure of 
gas mass, which we can obtain from this optically thin 
line by "counting molecules ". If we adopt abund ances 



Wilson fc Roodl[l994), and 



(e.g. |Strong et al 1988 Hunter et al. 1997Q, we 



would derive a molecular gas total mass of Mu 2 {R < 
500pc) = 1.6 x 10 9 M for the central kpc, which is 



of [ 12 CO/ 13 CO] = 60, (e.£ 

[ 12 CO/H 2 ] = 8.5 x 10- 5 ( IFrerking et al.||198"2^ ~and an 
excitation temperature, T ex = 20 K, which are typical of 
clouds in the nuclei of nearby galaxies, then N(H 2 ) col- 
umn densities can be estimated from 13 CO line strengths. 
The peak nuclear column densities derived from 13 CO are 
3.6 times lower than the corresponding values obtained 
from 12 CO using the Galactic Xqo- Our 13 CO mass 
suggests that Xg£ c ~ 0.6 x 10 20 cm~ 2 (K km s -1 ) -1 in 
IRAS 04296, similar to values found for other spiral nu- 
clei, LIRGs, and ULIRGs. 

We conclude that the nuclear H 2 mass (including He) 
for IRAS 04296 is M H2 (R < bOOpc) ~ 4.3 x 10 8 M 
(H# 2 ~ 600 M pc -2 ) for the inner kpc, with estimated 
uncertainties of ~50% due to unknown 13 CO excitation 
temperature and abundance. Caution should be used 
in comparing this value to masses determined for other 
galactic centers in the literature in which the Galactic 
Xco is used to determine the mass; those masses are 
probably systematically overestimated. 

Outside the nucleus, the CO intensity in the disk re- 
mains high in absolute terms. Normalized by distance 
squared, the CO intensity of IRAS 04296 averaged over 
its central 45 aperture is lar ger than any galaxy ob- 
served in the FCRAO survey (Young et al. |1995| ) out 
to its 30 Mpc distance. Over the central aremmute of 
the OVRO field of view a total of M gas (R < 30") = 
6 x 10 9 M of molecular gas is implied (for a Galactic 
conversion factor includin g He) . Even ado pting the lower 
M 82 conversion factor ( | Walter et al.|2002| ) for this rather 
normal appearing disk, 1KAS 04296 nas twice the molec- 
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ular mass of M 82, including all outflows and streamers. 
13 CO(1-0) is weakly detected towards the disk in IRAS 
04296 so further constraints on the validity of the con- 
version factor are testable here as well. At 10 radii, the 
values of N(H2) derived from 13 CO are within a factor of 
1.6 of what the standard conversion factor predicts, and 
along the arms at distances larger than 10 (particularly 
at the bar ends), the two column densities agree within 
the uncertainties. This suggests that Xqo m the disk of 
IRAS 04296 may be consistent with Galactic disk values. 
However gradients in excitation temperature, 13 CO opac- 
ity and isotopic abundances can affect the 13 CO mass by 
at least a factor of two, so this merits further study. 

The IR luminosity of IRAS 04296 is - 9.8 x 
10 10 



L Q (Sanders et al. 
central arcminute, IRAS 

-1 ^J2\-l 



2003). Averaged over the 
TJ32SB has L IR /L CO = 72 
L (K km s- 1 pc 2 )- 1 , L IR /M H2 = 18 L /M and 
MQWMffj) ~ 5, all typical of LIRGs (e.g. Scoville et al. 
1991 Sanders fc Mirabel|1 996). Given the rotation curve 
, an inclination ol 50° and the molecular gas masses 
derived above, molecular gas accounts for at least 30% 
of the dynamical mass within the central 500 pc radius 
(using the lower 13 CO mass) and 15% of the dynamical 
mass over the entire mapped region. 

4.3. Gas Stability Against in the Molecular Bar 

Is the molecular gas disk in IRAS 04296 so dense that it 
will collapse into new stars without further triggering, or 
is the disk stable? We can assess this using the Toomre 
Q parameter, 



shear/non-circular motions in the strong bar act to in- 
hibit star formation in the arms. The morphology of CO 
and the starburst's location at the inner terminus of the 
bar (Fig. [8]) argue for bar induced cloud-cloud collisions 
as inflowing molecular gas piles up in the nucleus. 

5. ENHANCED STAR FORMATION EFFICIENCY IN THE 
NUCLEUS 

From the nuclear mass and disk mass determinations, 
and the star formation rates (SFR), we can compare nu- 
clear and disk star formation efficiencies (SFE). We con- 
sider three representations of SFE: 1) the gas consump- 
tion timescale, tsf = (^gas/^SFR), 2) the ratio of the 
mass of formed stars to total mass available in stars + 
gas, r] = [M sta rs/(M gas + M sta rs)], and 3) L IR /M gas . 

We obtain the SFR for disk (0.5 kpc <R< 4.3 kpc, 
not inclination corrected) and nucleus (R < 500 pc) in 
different ways. We estimate the nuclear SFR from our 
free- free fluxes for the inner 5" region, as described in 



where n is th e epicyclic frequency (e.g. Binney 



Tremaine 1987), a is the gas velocity dispersion, H gas 
is the total gas surface density and a is a constant of or- 
der unity depending on the stru cture of the galactic disk 
(here assumed to be a = 1) ( |Safronov||1960| |Toomre 
|l964[ ). If Q~ l the gas disk is unstable. 

Figure flO] displays the observed values of T, gas (cal- 
culated both from 12 CO and 13 CO, neglecting the HI 
surface density contribution, which if conservatively 
estimated to be distributed uniformly over the disk 
contributes much less than 10% everywhere within 
R=17. 5), <7, k and Q as a function of galactocentric ra- 
dius. The data has been averaged in azimuth assuming a 
disk inclined at 50° to the line of sight. For 12 CO and the 
Galactic conversion factor, Q ~ 1 for the inner 5 and 
rises to 3 - 4 by the end of the bar. This dependence is 
driven primarily by the rapid drop in gas surface density 
just outside the nucleus together with the fact the veloc- 
ity dispersion of the gas remains fairly large (a ~ 20—40 
km s _1 ) along the entire molecular bar. At first glance 
this is consistent with gravitational instabilities driving 
the starburst. But note that this depends on the overes- 
timated H2 mass from the standard conversion factor. If 

we use T* gas derived from optically thin 13 CO emission, 20 times shorter than the tq^ 



£4.1 (which are also consistent with the 11.7/im LWS 
flux). We obtained 7Vj™ c c = 6 ± 1 x 10 53 s~\ equivalent 
to Uq-q = 2 x 10 10 L»0, or a total stellar mass asso- 
ciated with the young stars of M^ a c rs = 2 x 10 7 M , 
and an estimated SFR nnc = 12M yr- 1 . For the disk 
SFR, we cannot use our radio continuum maps, since 
they are not sensitive to extended emission. To estimate 
the SFR of the disk, we use the IRAS 12/im flux and 
subtract off the nuclear contribution using our 11.7/im 
LWS flux. About half of the IRAS 12/im flux originates 
from the inner 5" of IRAS 04296. The total disk SFR is 
only slightly higher than that of the nuclear starburst, at 
SFR di S fe ^ 13±5M yr- 1 , or 2-5 times that of the Milky 
Way. The total predicted luminosity in massive stars (if 
we assume the 12/im flux is excited by OB stars) pro- 
duces only half of the observed total infrared luminosity, 
which will also have contributions from both older pop- 
ulations of stars and regions of exclusively low mass star 
formation. Our SFR value for the disk may therefore 
be an underestimate of the true SFR, but only signifi- 
cantly so if the predominant mode of star formation is 
exclusively low mass, which we view as unlikely. 

To determine the nuclear SFE, we use the nuclear 
gas mass (from 13 CO including He) is M™s = 4.3 x 
10 8 M . The nuclear SFE can then be represented as 
(Lm/M gas ) nuc ~ (L OB /M gas ) nuc ~ 30-50, with the un- 
certainty due to the gas mass. This is high, but within 
the ra nge found for starburst galaxies ( [Young fc Scoville] 
1991). The surface density of star formation tor this re- 
gion is then £g^(R< 500 pc) = 15 ± 9 M yr" 1 kpc" 2 . 
The surface density of gas, which we assume is all molecu- 
lar, is S gas (R < 500pc) = 600 M pc~ 2 . This SFE point 
lies about an order of magnitude abov e the nominal re - 
lation for the Kennicutt- Schmidt law ( Ke nnicutt| [l998 ) , 
although within the scatter. The gas depletion timescale, 
T sW = 4 x 10 7 yrs, 50 times shorter than the disk value 
of t~sf = 2 Gyr foun d for the 23 local s piral galaxies of 
the THINGS survey (|Leroy et al]|2008[) and mo re than 



a rather different picture appears. In this case derived 
Q's would be effectively constant with radius over much 
of the bar and the disk would be stable everywhere. This 
suggests that the starburst is not a result of gravitational 
instabilities in a disk alone and possibly that the large 



"Gyr in M51 flKenni- 
cutt et a"T]|2007 Blanc et al. 2009) at similar gas surface 



densities. In terms ol the percentage of molecular mass 
converted to stars, n nuc c± 5 ± 1% over the 500 pc of our 
CO beam, and is almost certainly significantly higher on 
smaller GMC scales. This SFE is already somewhat high 
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for GMC scales of 50 pc or more, but is characteristic of 
SFEs seen in the Ga laxy on cluster (<10 pc) scales (e.g. 
Lada & Lada 2003| ). Therefore star formation appears 
to progress at an elevated, but not unusual efficiency in 
the nucleus. 

From the SFR and molecular gas mass estimated 
for the disk, we can compute the disk SFE. Using 
the Galactic value of the CO-H2 conversion factor ap- 
propriate for disks, we compute a disk-only molecu- 
lar gas mass of M^ k = 4.3 x 10 9 M , which is 10 
times the mass of the nuclear gas. The SFR over 
the entire disk, by contrast, is comparable to that of 
the nucleus, at SFR disk ~ 13 ± SA^yr" 1 ,. We find 
(L IR /M gas ) disk ~ (L OB /M gas ) d ™ k ~ 5, uncertain to 
about a factor of two largely due to gas mass. This is 
similar, wi thin the uncertainties, to the global Galactic 



Scoville|1991 , corrected to our value 



value of 8 ( |Young < 
of Xqq). The surface density of star formation is then 

^sfr ~ °- 3 ± °- 15 M © y r_1 kpc" 2 . The surface density 
of gas, which we assume is all molecular in this region, 
is E^f = 75 M Q pc- 2 . The disk SFE can also be rep- 
resented by TiSFR/^f'as = 3.0 x 10 -9 yr~ l , or a gas 
consumption timescale of r^ k ~ 3.3 x 10 8 yrs. The 
gas consumption timescale is a factor of six short er than 
the typical THINGS galaxy (Leroy et al. 2008), and a 
factor of three sh orter than the average value for M51 
( |Blanc et al.||2009| ). In terms of percentage of molecular 
gas converted into stars, rj dlsk = 0.5%. Within the uncer- 
tainties, this is consistent with the SFE of ~1% observed 



for Galactic disk star formation on GMC scales ( Lada & 
Lada||2003] |Ikeda et al.|[2007| |Higuchi et aT1|2009| | John- 
ston et al.||2009p. 



A low disk SFE is consistent with the model of a nu- 
clear starburst produced through secular evolution. If 
nuclear star formation is seen, it means that fresh gas 
must be present in the nucleus, else it would have formed 
stars long ago. In a morphologically undisturbed galaxy 
such as IRAS 04296, the nuclear gas must be a product 
of secular evolution. If the SFE of the disk were higher, 
then the gas would be used up before it arrived in the 
nucleus. Much of the molecular gas presently lies within 
1.7 kpc (12") of the center. At this radius the orbital 
timescale is ~ 10 8 yrs. If we estimate the velocity for the 
outermost gas in IRAS 04296 to dri ft inward via secular 
evolution to be a few km s _1 (§3.3 |Athanassoula||l992 



Kormendy fc Kennicu tt 2004), then the timescale for al 
of the CO gas in the disk of IRAS 04296 to drift into the 
nucleus is ~1 Gyr. The star formation timescale aver- 
aged over the disk is r^p k ~ 3 x 10 8 yrs. If the SFE were 
much larger than this, then there would be no build-up 
of nuclear gas, and we would not see strong nuclear star 
formation. 

We conclude that the SFE in the nucleus of IRAS 
04296 is an order of magnitude higher than the SFE in 
its disk. The value rj disk ~ 0.5%, for the disk of IRAS 
04296 is consistent with the Milky Way value of 77 ~1% 
on GMC scales. The value of r\ nuc > 5% found for the 
nucleus is consistent with SFE values on cluster (<10 pc) 
scales in the Galaxy. The corresponding star formation 
timescales are r\ 



—disk 
T SF 



SF 



3 x 10 7 yrs for the nucleus and 
3.3 x 10 8 yrs for the disk. While the disk perco- 
lates along at normal Galactic disk SFRs, the starburst 



mode "waits" until the gas has drifted to the nucleus. 

6. COMPARISON OF IRAS 04296 WITH OTHER NEARBY 
GALAXIES AND LIRGS 

IRAS 04296 is a LIRG of the local universe. We now 
place IRAS 04296 in the context of the other better- 
known members of that class: NGC 2146, NGC 7552, 
NGC 4418, NGC 1365 and NGC 1068. NGC 1068 and 
NGC 1365 contain active nuclei which contribute signifi- 
cant fractions of their total luminosity. NGC 7552 hosts 
a weak LINER. NGC 4418 is thought to be a Seyfert 
2 because it contains a very a comp act bright nucleus 
dEales et al.||1990| |Evans et al. |2003[ ) , but it is no more 



compact than the infrared continuum from IRAS 04296. 
NGC 4418 actually resembles a compact nuclear star- 
burst more than the composite galaxies with which it is 
usually classed. NGC 2146, along with the less luminous 
galaxies M82 and the Antennae, are thought, like IRAS 
04296, to be pure starbursts. 

In contrast to NGC 4418 and IRAS 04296, star for- 
mation is more spatially extended in the remainder of 
the LIRGs. NGC 7552 and NGC 1068 have kpc-scale 
starburst rings, well outside the nucleus. The starburst- 
dominated galaxies, NGC 2146, NGC 1365, M82, and 
the Antennae, have starbursts extended over many hun- 
dreds of parsecs. NGC 7552 and NGC 1365 are the only 
other starburst-dominated local LIRG that does not ap- 
pear overly disturbed; both have barred stellar appear- 
ances similar to IRAS 04296's. NGC 7552 may have 
experienced a similar dynamical history to IRAS 04296, 
namely slow, bar-mediated evolution, although in IRAS 
04296's case the starburst is concentrated in the nucleus, 
whereas NGC 7552's enhanced star formation occurs well 
before the gas reaches the nucleus. 

The local universe LIRGs are a surprisingly diverse 
class. Many of the differences are probably because they 
are at different evolutionary stages but others appear to 
be more fundamental. They may perhaps relate to the 
underlying causes and mechanisms of the starbursts. 



7. 



CONCLUSIONS: FIRST LOOK AT IRAS 04296, A 
BAR-DRIVEN STARBURST TURNING ON 



The little-studied nearby (29 Mpc) LIRG, IRAS 04296, 
is a barred and normal-appearing spiral galaxy, with an 
intense nuclear starburst and a vast reservoir of molec- 
ular gas. Mid-IR and radio continuum imaging show 
that IRAS 04296 currently has a bright nuclear starburst 
containing the equivalent of ~ 10 5 07 stars (N^y^ ~ 
6 ± 2 x 10 53 This ~150 pc region alone has the 

equivalent of more than half the star-forming luminos- 
ity in the entire 500 pc starburst region in M82. The 
nuclear starburst region accounts for ~50% of the to- 
tal current star formation in the galaxy as measured by 
11.7/im continuum flux, and contributes ~20% of the to- 
tal IRAS luminosity of IRAS 04296. The starburst lumi- 
nosity is concentrated in a double source that is <50 pc in 
size, located near the inner end of the northwestern bar 
arm. Radio and mid-IR images suggest that these are the 
HII regions associated with two massive young super star 
clusters containing a total of 2±1 x 10 T M of stars, and a 
total young stellar luminosity of Lob = 2±0.7x 10 10 L . 

The morphology and kinematics of the CO emission of 
IRAS 04296 show a starburst that is being fueled by slow 
gas inflow along a bar. We can estimate the history and 
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lifetime of the starburst from the gas content and inflow 
velocity. If the velocity of gas inflow to the nuclear star- 
burst is the maximum value of the average observed bar- 
induced peculiar velocity, ~20 km s _1 , a net radial inflow 
of gas along the arm width (~750 pc at 10 ) of 9 M yr _1 
results (for a Galactic conversion factor relevant to the 
arms). For more expected net inflow rates of a few per- 



cent of the orbital velocity ( Athanassoula 1992), M in f < 



1 M yr _1 . The star formation rate in IRAS 04296 from 
the thermal radio continuum is consistent with a current, 
nuclear star formation rate of SFR nnc ~ 12 ± 6 M yr _1 
and as much as SFR £o£ - 25 ± 10 M yr" 1 total for the 
galaxy as a whole. The nuclear star formation is not in 
equilibrium with the net inflow rate and hence the star- 
burst is currently consuming its reservoir of nuclear fuel 
faster than it is being replenished. 

Stellar mass buildup also dictates that the current rate 
of nuclear star formation is not sustainable over the long 
term, that it is a true starburst. Subtracting the gas 
mass from the dynamical mass leaves a small nuclear 
stellar mass of <1.2 x 10 9 M . Nuclear star formation 
also could not have been occurring at the current ~12 
M yr _1 rate for more than the past ~ 10 8 yrs or the 
stellar mass would exceed the observed dynamical mass. 

We also find that the star formation efficiency is a func- 
tion of position in the galaxy. We derive a starburst 



timescale of r^p c 



3.7 x 10 T yrs for the nucleus (R 



< 500 pc), or r] nuc =b%] the star formation timescale 



is ten times longer in the disk, and the SFE is only 
77 nuc =~0.5%. The relatively low SFE observed in the 
disk, where the gas consumption timescale is ~5 orbital 
times, is necessary so that gas drifting into the nucleus 
along the bar orbits would not be consumed before it 
reached the center. Any higher value for the disk SFE 
would be inconsistent with a nuclear starburst fueled by 
gas drifting inward via bar inflow. 

Taken together these facts imply that we are witnessing 
an early stage of starburst / bar inflow induced secular 
evolution in the history of IRAS 04296. Otherwise the 
bulge would be much larger and the disk depleted of its 
gas. 
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TABLE 1 
IRAS 04296+2923 Basic Data 



Characteristic Value Reference 



Dynamical Center a a(2000) = 04^32 m 48 s .65 ± l" 1 

(kinematic) 5(2000) = +29°29 / 57. // 45 ± l" 

2/rni peak (2MASS) a(2000) = 04^32 m 48 s .60 db 0."3 1 

5(2000) = +29°29 / 57. // 49 ± 0." 3 

Y isr a 2086 kms- 1 1 

Adopted Distance 29 Mpc 1 

Position Angle a 252° 1 

Inclination 50° 1 

Rmacc(J band) 8.1 kpc 1 

Rmax (rotation curve) a 2.15 kpc 1 

Vmax a 190 km s _1 1 

n a 1.25 1 

ba r h "85° 1 

Q bar h 43 km s- 1 kpc" 1 1 

M# 2 (< 500 pc) d 4.3 x 10 8 M 1 

E^ 2 (< 500 pc) d 550 M pc" 2 1 

M dyn (< 500 pc) a 1.6 x 10 9 M 1 

M# 2 (<4.3 kpc) c 5.9 x 10 9 M 1 

M dyn (< 30 kpc) a 3.0 x 10 10 M Q 1 

Mgf 1.38 x 10 9 M 3 

IRAS 12/i, 25/z, 60/x, 100/x 1-39, 5.90, 42.1, 48.3 Jy 2 

Ljfl 9.8 x 10 10 L 2 



References. — (1) This paper; (2) Sanders et al. (2003); (3) Chamaraux 
et al. (1995) 

a Based on the best fitting Brandt rotation curve (§??). 

b Based on the best fitting bar model (Figure (9|. 

c Assuming the stan dard CO conversion factor^ ^4.2| . 
d From I3co(l-0) (j}4~2). 



TABLE 2 
Infrared and Radio Fluxes 



A r.m.s. Beam (Robust=0) 6 m ax Peak Flux Total Mapped a ' b Cut Flux b ' c Cut Flux b ' c 

(u,v)min = 18kA (u,v)mi n = 50kA 

mjy/bm " x" , pa° " mJy/bm mjy mjy mJy 



20 cm 


0.11 


1.60x1.35, 82.8 


38 


32 


130 ± 10 


71 ± 5 




6 cm 


0.04 


0.48x0.39, -83.3 


10 


3.8 


59 ± 5 


28 ± 3 




3.6 cm 


0.034 


0.25x0.21, 88.2 


7 


1.5 


25 ± 5 


14 ± 2 


11 ±2 


2 cm 


0.085 


0.18x0.138, -53.9 


4 


1.0 


10 ± 1 




10 ± 2 


1.3 cm 


0.068 


0.1x0.09, -80.39 


2 


0.5 


7 ± 1.5 




7 ± 1.5 


1.3 cm BC array 


0.16 


0.91x0.30, 76 


20 


3.7 


23 ± 5 


18 ± 3 




2.7 mm 


0.75 


4.6x3.8, 8 


52 


9.9 


11 ± 2 






11.7 /im 


5.1 cts/s 


0.4 


n.a. 


71±10 


680 ± 100 







a Based on "uncut" maps with full native A configuration (u, v) coverage. 

b Quoted uncertainties in the fluxes are based on 3 times the rms times \/]V times the number of beams corresponding to the source size; 
they do not include resolved-out flux, or the estimated 5% (10% at 2 and 1.3 cm) absolute flux scale uncertainty. Large beam fluxes do not 
exist for wavelengths other than 20 cm, so we cannot estimate amounts of resolved-out flux except at 20 cm, where it is negligible. Fluxes 
are measured within X/Bmin for each map. 

c Based on "cut" maps restricted in (n, v) coverage. The 18kA fluxes correspond to the maps of Fig. [3] and the 50kA fluxes to those of 
Figs. [5] and [6] See p3] 
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Fig. 1. — Left) The J band (1.25/im) image of IRAS 04296 taken on the Palomar 5-meter. Right) Palomar continuum-subtracted Pa 
/3 image of IRAS 04296. The many point sources in the image are foreground Galactic stars of the Taurus region. Residuals from the 
continuum subtraction are seen in Pa (3; in particular a large residual contribution remains from the bright foreground star to the northeast 
of the galaxy. These residuals are likely due to slight mismatches in the point spread function of the (partially) saturated J band stellar 
images. The scale bar is 20 /; (2.8 kpc) in each plane, and the orientation is north-up, east-left. 




Fig. 2. — The Keck LWS image of IRAS 04296 at 11.7/xm. The scale in counts per pixel is above the figure and the conversion is 
4.1 x 10 -5 Jy/count for a 360 second integration. North is up and east is left. 
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Fig. 3. — Maps of the radio continuum, emphasizing extended emission, (top left): 20 cm. (top right) 6 cm. (bottom left): 3.6 cm. 
(bottom right): 1.3 cm with the BC array. Contour intervals are ±2 n / 2 (3cr) of the following increments: 0.25 mJy, 0.15 mJy, 0.12 mjy, 
and 0.25 mJy/bm for 20, 6, 3.6 and 1.3 cm, respectively. The dynamical center is displayed in the 20 cm and 6 cm maps. 
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Fig. 4. — The spectral index between 1.3 and 20 cm in greyscale, overplotted with representative contours of the 20 cm emission. The 
spectral index is not uniform over the whole frequency range covered, but maps at other frequencies appear similar to this. 
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Fig. 5. — Matching 3.6, 2, and 1.3 cm continuum images of IRAS 04296 made from "cut" data as described in the text, with consistent 
(u,v) coverage and a beamsize of C//15 x / /14, p. a. -66°. The bright central source, elongated east-west with a size of ~0'/5 has a flat 
spectum, containing flux ~ 3 mJy. The central arcsecond contains ~6 mJy at 1.3 cm. Contours are ±2 0-52 x 0.5 mJy/beam (~2.5-3<r at 2 
and 1.3 cm). The 3.6 cm maps are most affected by the cut of short (it, v) spacings. a) (top left) 1.3 cm. Peak flux 1.5 mJy/beam, total 
flux 7 mJy. b) (top right) 2 cm. Peak flux 1.1 mJy/beam, total flux 10 mJy. c) (bottom left) 3.6 cm. Peak flux 1.6 mJy/beam, total flux 
10-11 mJy. d) (bottom right) 3.6 cm map as in c) except plotted with a lowest contour of 0.2 mJy/beam. Cross marks the location of the 
dynamical center from the CO maps. 
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Fig. 6. — Matching 2 and 1.3 cm continuum images of IRAS 04296. Contours are ±2°- 5i x 0.35 mJy/beam (2.5a). Beam for both 
maps is 0'/15 x 0"14, p. a. -66°. Left) 2 cm image, with (u,v) coverage cut to match the 1.3 cm image as described in text. Peak flux 0.9 
mJy/beam. Right) 1.3 cm map convolved to the same beam as the 2 cm image in a). Peak flux 0.8 mJy/beam. 




Fig. 7. — Far left) The 2MASS K s image of IRAS 04296 over the same region covered by the molecular data. Contours are logarithmic 
in arbitrary units. Middle Left) The naturally weighted 12 CO(1-0) intensity map. Contours are 13.9 K km s _1 x 1,2,3, ...20, 30,40, 50 for a 
beamsize of 4. 6 x 3. 6; -14°. Middle right) The 13 CO(1-0) integrated intensity map of IRAS 04296 displayed at the same resolution as 
12 CO(1-0). Contours are in steps of 10.3 K km s — 1 . Far right) The 2.7 mm continuum map generated from the line free channels of the 

wide band dataset. Contours are in steps of 2.0 mJy beam -1 (2a) for a beamsize of 4. 6 x 3. 8; —15°. The crosses in the middle planes 
mark the location and extent of the 6 cm radio continuum emission. 
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Fig. 8. — 6 cm radio continuum (bold lines) overlaid on the central region of the robustly weighted CO (1-0) image (thin lines). CO (1-0) 
contours are in steps of 4 Jy bm _1 km s _1 , while the 6 cm contours are 0.2 mJy bm _1 x -3,-1,1,3,5.... The 6 cm beam is in the lower left 
and the CO(1-0) beam (2." 6 x 2."l; 21.6°) in the lower right. 
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Fig. 9. — An epicyclic weak-bar analytical model for IRAS 04296. Upper left) The CO(1-0) map (see Figure [Tt. Upper center) The 
corresponding model for the morphology of the gas. Regions of high dot concentration reflect the expected locations of gas populated orbits. 
Upper right) The same model de-inclined to show the face-on morphology of the bar. Lower left) The Posit ion- Velocity (PV) diagram 
taken along the major axis (see Table Q of the galaxy, uncorrected for inclination. The zero velocity coordinate in the figure corresponds 

to 2079 km s" 1 and the zero position coordinate corresponds to a(J2000) = 04 h 32 m 48. s 56 ; 5(J2000) = 29° 29' 57." 6. Crosses mark 
the fitted azimuthally averaged rotation curve. Uncertainties are given by the size of the symbol except for the points at the largest radii, 
which have significantly higher uncertainties. Lower right) The PV diagram for the model, generated in a consistent fashion for comparison 
with the observed PV diagram. 
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Fig. 10. — The surface densities and Toomre Q parameters for IRAS 04296. The dot-dashed lines are the azimuthally average molecular 
surface densities, S# 2 , as a function for galactocentric radius. Black lines correspond to S pas calculated using CO(1-0) and t he G alactic 
CO-to-H2 conversion factor, while gray lines correspond to ^ ga s calculated from 13 CO(1-0) assuming it is optically thin (j ]3.3[ ), both 
inclination corrected. The velocity disperson a (dashed line) is taken from the line width of the CO (1-0) transition, and n (dotted line) 
from the fitted rotation curve (§??). The resulting Toomre Q parameter (j |4.3) is displayed as the thick solid line (black for values from 
CO determined surface densities and gray for 13 CO gas surface density determinations). Units and scaling for the vertical axis are given 
in the figure legend. 



